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0\ ■ ABSTRACT 

^ I We present the geometrical optics for refraction of a distant background radio 

' source by an interstellar plasma lens, with specific application to a lens with a 

^ ■ Gaussian profile of free electron column density. The refractive properties of the lens 

are specified completely by a dimensionless parameter, a, which is a function of the 



lazio@rsd.nrl.navy.mil 



wavelength of observation, the free electron column density through the lens, the 
lens-observer distance, and the diameter of the lens transverse to the line of sight. 
5^ I A lens passing between the observer and a background source, due to the relative 

motions of the observer, lens, and source, produces modulations in the light curve of 
the background source. Because plasma lenses are diverging, the light curve displays 
a minimum in the background source's flux density, formed when the lens is on-axis, 
surrounded by enhancements above the nominal (unlensed) flux density. The exact 
form of the light curve depends only upon the parameter a and the relative angular 
sizes of the source and lens as seen by the observer. Other effects due to lensing 
include the formation of caustic surfaces, upon which the apparent brightness of the 
background source becomes very large; the possible creation of multiple images of the 
background source; and angular position wander of the background source. If caustics 
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are formed, the separation of the outer caustics can be used to constrain a, while 
the separation of the inner caustics can constrain the size of the lens. We apply our 
analysis to two sources which have undergone extreme scattering events: 0954+654, a 
source for which we can identify multiple caustics in its light curve, and 1741—038, 
for which polarization observations were obtained during and after the scattering 
event. We find general agreement between modelled and observed light curves at 
2.25 GHz, but poor agreement at 8.1 GHz. The discrepancies between the modelled 
and observed light curves may result from some combination of substructure within 
the lens, an anisotropic lens shape, a lens which only grazes the source rather than 
passing completely over it, or unresolved substructure within the extragalactic sources. 
Our analysis also allows us to place constraints on the physical characteristics of the 
lens. The inferred properties of the lens responsible for the scattering event toward 
0954+654 (1741-038) are that it was 0.38 AU (0.065 AU) in diameter, with a peak 
column density of 0.24 pc cm^^ (10"'* pc cm~^), an electron density within the lens 
of 10^ cm-3 (300 cm-3), and a mass of 6.5 x 10"^^ Mq (10"^^ Mq). The angular 
position wander caused by the lens was 250 mas (0.4 mas) at 2.25 GHz. In the case of 
1741—038, we can place an upper limit of only 100 mG on the magnetic field within 
the lens. 



Subject headings: astrometry — ISM:general — radio continuumdSM 



1. Introduction 

Strong refraction by interstellar electron density irregularities has been invoked to explain 
several radio astronomical phenomena. Fiedler et al. (1987, 1994a) and Clegg, Fey, & Fiedler (1996) 
reported a number of extreme scattering events (ESEs) — dramatic changes in the 2.25 GHz fiux 
density — in the light curves of several extragalactic radio sources observed during the course of a 
2.25 and 8.1 GHz monitoring program. The changes typically lasted several weeks or months and 
were occasionally accompanied by strong variations at 8.1 GHz. Because of the simultaneity of 
the events at different wavelengths, the time scales of the events, and speed-of-light travel time 
arguments, Fiedler et al. (1987) and Romani, Blandford, & Cordes (1987) concluded that the flux 
density variations were likely due to strong scattering by the Galactic interstellar medium (ISM). 
A similar dramatic change in the flux density, and times of arrival, has since been observed during 
a timing program of the pulsar PSR B1937+21 QCognard et al. 1993|) . 



Interstellar refraction may also be responsible for episodes of multiple imaging and fringing 



observed in various pulsars' dynamic spectra (Hewish, Wolszczan, & Graham 1985; Wolzszcan &: 



Cordes 1987| ; Gupta, Rickett, & Lyne 1988; Clegg, Fiedler, & Cordes 1993). Dynamic spectra 
record the flux density of the pulsar as a function of time and frequency. Because a pulsar is an 
angularly compact coherent emitter, strong interstellar refraction may produce multiple images of 
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the pulsar. The multiple pulsar images can then interfere with one another, producing fringes in 
the time and frequency domains which are visible in the dynamic spectra as alternating bright 
and dark bands modulating the diffractive scintillation pattern of the pulsar. 

The characteristics of the medium responsible for strong interstellar refraction are not clear. 
For example, it is unknown if the refraction responsible for such phenomena as ESEs and fringing 
in pulsar dynamic spectra is the result of a localized large-amplitude electron density enhancement 
(a "plasma lens") located somewhere between the observer and source, or whether the refraction 
is the result of an extended turbulent medium with a large number of electron density fluctuations 
distributed along the line of sight. There are arguments that favor both models: Clegg, Chernoff, 
& Cordes (1988) showed that expanding shock fronts can naturally create discrete plasma lenses, 
while the general ISM is widely believed to have an extended, turbulent distribution of electron 
density fluctuations ( Rickett 19901) . 



To date, refraction phenomena which have been interpreted in terms of the discrete lens model 



have been analyzed using back-of-the-envelope calculations for the effects of the lens (e.g., Clegg 



et al. 199^ ). These calculations have been useful for deriving approximate characteristics of the 



background radio source and the lens itself. The purpose of this paper is to venture beyond the 
back-of-the-envelope calculations and to present a detailed, quantitative description of the optics 
properties of an interstellar plasma lens. The motivation is to present the optics of interstellar 
plasma lenses in a convenient format that allows the application to a wide variety of problems. 

In this paper, the analysis is limited to the refractive properties of a one-dimensional plasma 
lens. Results are applied to the specific case of a lens with a Gaussian distribution of free electron 
density across the lens plane. The Gaussian lens is a convenient form for describing a plasma 
irregularity confined to a characteristic size, and whose influence on the background wave front is 
negligible outside of this characteristic size. 

The plan of this paper is the following. In §^ we present a general description of the 
geometrical optics of a plasma lens, in §|3| we specialize the application to the case of a Gaussian 
lens form, in we discuss the light curves expected from a Gaussian plasma lens passing in front 
of both point and extended sources, in §|5| we discuss the formation of caustics and their utility 
for extracting information about the lens from the light curve, in §^ we discuss the expected 
angular position wander caused by a Gaussian lens, and in ^ we apply our model to the extreme 
scattering events seen toward the sources 0954-1-658 and 1741—038. We summarize our conclusions 
in the last section. 



2. Geometrical Optics of a Plasma Lens 

The radio frequency index of refraction in an unmagnetized plasma is 
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= l-A^rene/vr, (1) 

where = Airnee'^ / nie is the square of the electron plasma frequency; e, me, = e^/(meC^), 
and rig are respectively the electron charge, mass, classical radius, and local number density, 
(wc = 56s^^y^ for rig in units of cm~^); lo = 2ttc/X is the angular frequency of the wave; A is its 
free-space wavelength; and c is the speed of light. 

The phase velocity through the lens is = c/rij-, which is always greater than c. Compared 
to the time taken by a wave in vacuum, the phase will traverse a lens of size I in a time 
T = l{l/c — less. In the limit cjg <^ uj, which is always true in the case of radioastronomical 

observations, the emergent phase from the lens will be advanced by an amount 

(f> = ojT = ArgiVe (2) 

compared to a wave in vacuum, where Ne = n^l is the column density through the lens. Because 
the phase is advanced, plasma lenses are diverging. 

Allow a one-dimensional plasma lens to have a profile of column density Ni.{x); implicit in 
this discussion of a plasma "lens" is the notion that Nf>{x) varies across the lens in a regular 
fashion. A plane wave from an infinitely distant point source is incident on the lens with phase (pQ. 
The phase as a function of x on the emergent side of the screen will be 

(/>(x) = 0o + ArgiVe(x). (3) 

In the geometrical optics limit, rays of energy flux propagate normal to the surfaces of constant 
phase. As A'g varies across the lens, there will be a slope to the emerging constant phase surfaces. 
The rays are therefore refracted by the lens through an angle 9^- ( Born fc Wolf 1980| , Chapter III): 



The observer is a distance D from the lens, Fig. ^a. The coordinate in the observer's plane is 
given by x' . A ray emerging from the plasma lens at coordinate x is refracted by an angle 6r{x) 
and strikes the observer's axis at a position 

x' = x- 9,{x)D, (5) 

for 9^{x) <^ 1, which will be the case throughout this paper. Multiple solutions x = x^ (with 
A: = 1, 2, . . . , n) to equation (|5|) can exist for a given x' . In such cases, an observer at x' sees n 
images of the background source appearing to come from the directions towards x\^ X2, • . . , ^n* 
The brightness of each image will be enhanced (or reduced) by some factor due to refractive 
focusing (or defocusing). The factor is determined by the spreading or focusing of rays due to 
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the gradient in refraction angle across the lens: 



Gh — -r-. 



dx 
d^' 



l-D 



de. 



dx 



(6) 



where solution to equation (^), is ultimately dependent on x\ D, A, and the function 

Nc{x), so that Gk = Gk[x',D,X,Nc{x)]. 

The same analysis can be modified slightly to account for an angularly extended background 
source. In this case, a source of angular extent is considered to be the superposition of 
point sources distributed over a range of angles — ^s/2 < 9 < +9s/2. For an infinitely distant 
source the distribution of rays incident upon the lens screen has a range of incidence angles 
— < < +9s/2. The intensity of the rays coming from each incidence angle Oi is weighted by 
the function B{9\) to account for the brightness distribution across the source. 



For rays incident upon the lens screen with incidence angle 9\, we obtain 



X 



x - [Orix) + e;]D. 



(7) 



The term 6iD in equation (^) accounts for the skewness of the incident ray with respect to the 
lens axis, Fig. The gain factor Gk is now also dependent on 9i: Gk = Gk[x\ D, 9i, A, Nc{x)]. 

The total observed intensity of the background source is obtained by integrating the product 
of the gain factor G^ and the source brightness distribution B(9) over all 9, and summing the 
result of this integration over the n images of the background source: 



/■ + 00 

I[x', D, A, N,{x)] = / B{9)Gk[x', D, A, 9, N,{x)] d9. 



(8) 



3. The Gaussian Plasma Lens 

The basis of our analysis of lensing is studying the behavior of I[x\ D, A, Nc{x)] in the case 
where Ne{x) follows a Gaussian distribution in x. The Gaussian distribution is a convenient 
function for describing a lens that is localized to some region of characteristic size a and whose 
refractive strength declines quickly outside of this region. In this section we develop the refractive 
properties of the Gaussian lens in the limit of geometrical optics. We will show that these 
properties can be parameterized in terms of a single dimensionless quantity that is constructed 
from the physical parameters A, a, D, and Nc(x). 

The column density profile takes the form 

N^{x) = Noexp\-{x/af]. (9) 
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Following the discussion in §|2|, the phase advance due to such a lens is 

(f){x) = XrcN^exp —{x/af' . 
The refraction angle, equation becomes 

dr{x) = —{X^reNo/Tra'^)xex.p —{x/af' 
and the ray path, equation (|^), is 

X =x\l + {X^r^NoD/ira^) exp [-(x/a)^] } - OiD. 



(10) 



(11) 



(12) 



Considerable simplification takes place if the following dimensionless variables are defined: 



u 



7 



x/a; 
x' /a; 
a/D- 

u' + /3s. 



(13a) 
(13b) 
(13c) 
(13d) 
(13e) 



The angle 6^ is roughly the angular width of the lens as seen by the observer, while /3s is the 
incidence angle of the background rays upon the lens screen in units of 9^. 

In the above equations the physical parameters A, a, D, and Nq occur consistently in the 
combination a. 



a = X^rcNoD/TTO^, 
\ a IT 



(14) 



We have written a in this second form to emphasize the essential physics. The phase advance 
through the lens is XtcNq/tt. The Fresnel scale is V XD. Thus, the properties of the lens are 
determined by the square of the ratio of the Fresnel scale to the lens size and the phase advance 
through the lens. The larger the parameter a, the greater are the observable effects due to the 
lens. Hence a weak lens, AreA^o/vr <C 1, can produce large observable effects if the Fresnel scale 
is sufficiently larger than the lens size. Similarly, a strong lens, XreNo/ir <C 1, need not produce 
large observable effects if the Fresnel scale is small relative to the lens size. Numerically, 



a 



3.6 



A 



No 



D 



-2 



1 cm 



1 cm 3 pc / \ 1 kpc / \ 1 AU 



(15) 



Upon substitution of the above dimensionless variables, the following expressions are derived 
for the refractive properties of a Gaussian lens: 
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refraction angle 
ray path 
gain factor 
total intensity 



6»,(u)/6'£ = -auexp(-u2); (16) 

u[l + a exp(-n2)] -7 = 0; (17) 

Gk = [l + {l-2ul)aeM-4)r^; (18) 

" r+oo 

I{u',a) = Y, B{P,)Gk{u',a,Ps)dps- (19) 

1 1 J —oo 



k=l- 



4. Characteristic Light Curve Produced by a Gaussian Lens 

When transverse motion between source, lens, and observer is assumed, e.g., motion 
of the observer along the u' axis, I{u',a) translates into a light curve I{u' ,a;t,v) since 
u'{t) = u'{t = 0) + vt where t is time and v is the relative transverse velocity. We will show 
that interpretation of radio light curves in terms of plasma lenses that can be approximated as 
Gaussian in profile can lead to inferences regarding the physical properties of the lens. 



4.1. General Description 

We will present numerical results in following sections. Here we will first develop a physical 
understanding for the basic characteristics of refraction by a Gaussian plasma lens, Fig. ||. Plane 
waves from an infinitely distant source are incident on the lens screen. The plane waves are 
indicated by straight dotted lines in the figure, and the lens is represented by a plot showing the 
electron column density as a function of coordinate u along the lens plane. This representation is 
purely schematic, as the lens is assumed to have a negligible but uniform width along the line of 
sight. Upon emergence, the constant phase surfaces are distorted into contours that mimic the 
function Ne{u), as represented by the inverse Gaussian shaped dotted lines. The lines are inverse 
Gaussian shaped because the phase velocity is greatest through the center of the lens, where A'e 
reaches its maximum value of Nq. 

In the limit of geometrical optics, the rays of energy fiux travel perpendicular to the constant 
phase surfaces. The direction of travel of the rays is indicated schematically by the small arrows 
in the figure. Extending this concept, we have drawn the path of approximately 100 rays as they 
travel perpendicular to the constant phase surface after emergence from the lens. The ray path 
is shown for a total distance D along the lens axis. An observer located close to the lens plane 
but far from the lens axis (e.g., the upper left and upper right regions of the ray trace) sees an 
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unchanged source: There is no change in the observed flux density (no change in the number 
density of the rays) or in the source's position (the direction of arrival of the rays). 

Somewhat farther from the lens and slightly off-axis, in the "Focusing Regions" marked in 
Fig. |2|, the rays begin to converge and their number density increases. An observer located in this 
area of the ray trace would see a source of enhanced brightness somewhat displaced from its "true" 
position, as evident from the increased number density and skewness of the rays, respectively. At 
the same distance from the lens, but closer to the lens axis, the rays are spread apart by the lens. 
An observer in this region would see a source of decreased brightness due to the lower number 
density of rays. The source would be offset from its true position by an amount that increases 
with the observer's distance from the lens axis. Directly on the lens axis, the number density of 
rays reaches a minimum, but the source appears in its correct position. 

At a larger distance from the lens plane, the skewed rays that make up the focusing regions 
cross one another. Inside this region where the rays have crossed, more than one bundle of rays 
would be incident on an observer. Each bundle would appear to come from a different direction. 
The observer therefore sees more than one image of the background source. In the case of a 
Gaussian lens, the greatest number of images seen is three. Since each ray bundle arrives from a 
different part of the lens, the spreading of the rays in each bundle, and therefore the brightness of 
the particular image represented by each bundle, may differ. 

Outside the focusing region, only a single image is seen. At large distances from the lens axis, 
the image has its nominal brightness, as evident from the uniform density and straight path of the 
rays in this region. Close to the lens axis the image has a reduced brightness as evident from the 
spreading of rays in this region. 

Separating the single image regions from the multiple image regions are the caustic surfaces. 
On these surfaces, an observer sees multiple images merge into a single image, and the local 
number density of rays from that image grows without bound in the limit of geometrical optics. 
In the case of a Gaussian lens, two of the three images merge on the caustic surface. In reality, the 
brightness of the image is limited by diffraction effects which are not taken into account here. 

Now assume that an observer, due to relative transverse motion of the lens, travels at a 
constant (large) distance D from the lens plane. The transverse coordinate along the observer's 
path is given by u' , which is a function of time, t. As the observer encounters different ray 
densities, the brightness of the background source (and its apparent position) will change with 
time. At the bottom of Fig. ^ is a light curve constructed by plotting the number density of 
rays encountered as the observer travels across the bottom of the ray trace. The intensity / is 
normalized to unity in the absence of the lens. In the multiple imaging regions, the observed 
brightness is taken as the sum of the number densities of the ray bundles that reach the observer, 
in accordance with equation (^). 

When the observer crosses a caustic surface, the number density of rays tends towards infinity 
(for geometrical optics), and therefore the observed brightness of the source becomes very large. 
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In our refraction simulations, this increase in intensity is represented by tall spikes in the light 
curve. The height of the spikes is generally not properly represented in refraction simulations, but 
qualitatively their existence is correct. 



4.2. Numerical Results 

4-2.1. Point Source 

We start by considering the refraction of a background point source by a Gaussian lens. In 
this case /3s = 0, and we take Iq to be the intensity of the background source in the absence of 
lensing. The observed intensity, equation (|^), reduces to 

n 

Iiu',a) = IoJ2Gk{u',a). (20) 

k=l 

For a given value of a, the function I{u',a) is determined in practice by numerically computing 
the root(s) Uk to equation (^) (which is a transcendental function) for a particular u' , then 
computing the gain factor(s) of equation (18) and performing the sum in equation (20). This 



is done for a range of u' to construct the light curve. 

As an initial example, consider a Gaussian plasma lens with peak electron column density 
Nq = 0.1 pc cm~^ and width a = 2 AU. The observer is located a distance D = 1 kpc away 
and is observing at a wavelength A = 20 cm. For this particular combination of physical 
parameters, a = 36, equation ([l5|); other combinations can also produce similar values of a. 
Equations (|l^)-(|l9|) have been solved over the range —20 < u' < 20 for this value of a, and 

Gk for each u' is plotted in Fig. ^. The plot represents the enhancement (or decrease) in the 
brightness of the background point source over its value in the absence of a lens. The features in 
the light curve follow the general discussion in the previous section. On the lens axis, u' = 0, there 
is one image and the intensity is at a minimum, with a gain factor of Gi = 0.027 in accordance 
with 1/(1 + a) as predicted by equation (|l8|). 



Other values of a will give different light curves. Larger a gives rise to greater transverse 
displacement of the rays. This can occur by (1) increasing the distance D between observer and 
lens; (2) increasing the maximum column density A^o or decreasing the characteristic size a so that 
the transverse gradient of the column density, and therefore the refraction angle, becomes greater; 
or (3) observing at a longer wavelength, where the index of refraction rir becomes smaller (greater 
deviation from 1). 

Light curves are shown for various values of a in Fig. ^. Consider first the cases of weak 
refraction, a = 0.1 and a = 1. The transverse displacement of the rays is small and no ray crossing 
occurs. Modulations in the light curve are due solely to focusing and defocusing of the ray bundles 
since there are no caustics. When a is increased to 5, the onset of ray crossing is evident from the 
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existence of caustic spikes. However, the observer is close to the distance D where ray crossing 
just begins to occur. This is evident from the close proximity between the inner and outer pairs of 
caustics, compared to the width of the defocusing region near the lens axis (cf. Fig. ^). We will 
say more in §^ on the distances between caustics and what that tells us about the lens. Moderate 
refraction is shown by the example a = 10. The pairs of caustics have moved farther apart and 
the minimum in the light curve on the lens axis becomes deeper: 1/11 of its nominal value. The 
light curve for a = 40 exhibits strong refraction. The outer pair of caustics forms very far from 
the lens axis, near the limits of the plot, while the inner pair is not as strongly affected. This fact 
will be exploited in §|5| to estimate the size of the lens based on observed light curves. On-axis, the 
source is only 1/41 of its nominal brightness. 

4-2.2. Extended Source 

For an extended source the effect of the lens is determined not only by a, but also by the 
relative angular extent of the source compared to the angular extent (~ a/D) of the lens. Since 
we will generally invoke a Gaussian brightness distribution across the source, the parameter /3s is 
taken as the ratio of the FWHM angular width of the source to the observed angular FWHM of 
the lens. 

The numerical simulations are shown in Fig. ^. Using an a = 25 lens we show the light curve 
from a point background source, /Sg = 0, as well as those for /3s = 0.1, 0.25, 0.5, 1.0, and 5.0. It 
is evident that an extended source smoothes out the caustic spikes and the transition into the 
defocusing region near the lens axis. However, the brightness minimum at the center of the lens 
remains substantial until the source is larger than the lens. For instance, when /3s = 1 the on-axis 
gain is 0.039, only 2% larger than the gain 1/(1 + a) predicted by a point source. When the 
source is much bigger than the lens the result of refraction is to produce a ripple in the light curve 
rather than strong modulations. Rays from the outer edges of the source are refracted into the 
defocusing region so the minimum intensity is not as deep as for a point source. 

Note that the height of the ripple and the true effects on the caustic spikes of an extended 
source can not be quantitatively determined without a simulation that also includes diffraction 
effects. We have performed limited diffraction simulations and can say that the refractive light 
curves are correct in a qualitative sense. 

5. The Caustics 

The caustics are a powerful probe of the physical properties of the lens. In our refractive 
simulations the important observable is the distance between the pairs of caustic spikes. Diffractive 
fringes form in the vicinity of the caustics due to interference between the multiple images. 
Though it is beyond the scope of this paper, the fringes can serve as a useful probe of both the 
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lens and the background source (e.g., Wolzszcan &: Cordes 1981^ ) . 



As evident from Fig. the distance between the inner two caustics measured in units of 
the characteristic size a of the lens is a relatively weak function of a, while the distance between 
the outer two caustics depends more strongly on a. We have solved for the locations where the 
caustics intersect the w'-axis, as a function of a. We define the dimensionless quantities Au( and 

as the distances between, respectively, the inner and outer pairs of caustics, in dimensionless 
units of u' . These quantities are plotted as functions of a in Fig. where we have also plotted 
(Auq — Aii-)/2, the distance between the outer and inner caustics on each side of the lens axis. It 
is evident from the figure that no caustics form if log^g a < logio Omin — 0.35 (amin — 2.25). In this 
case the observer is closer to the lens than the point at which ray crossing occurs. For a = amin, 
caustics have formed but the inner and outer caustics are merged together. As a increases, the 
outer pairs spread apart while AnJ remains relatively fixed. 

Although these separations are derived from solutions to transcendental equations, it is 
possible to find an analytic expression for An^ as a function of a, the derivation of which we 
present in the next section. For Au- we have performed a linear least-squares fit in the log-log 
domain to find an approximate functional dependence on a. The separations are 

An; ^ 4.30q°-°^^ (21) 
A< = ^/2(l + ae-i/2) . (22) 

The dimensionless quantity An| is Ax|/a, where AxJ is the distance between the inner 
two caustics. The characteristic size a of the lens can therefore be estimated as a ~ Ax;/4.3. 
Additionally, the value of a can be estimated by measuring the relative separations between the 
inner and outer pair of caustics: 

\ 1.06 



where Ax^ is the distance between the outer pair of caustics. 

We emphasize that the distance between the inner pair of caustics is the only reliable estimate 
for the size of the lens. Using the value Ax^ determined through observations of the total duration 
of the lensing event to judge the size of the lens is not appropriate, cf. Fig. |6[ It can deviate by a 
factor of order a from the true characteristic size a of the lens. Since physically plausible values 
of a can exceed 10^, a characteristic lens size derived from the total duration of the lensing event 
will overestimate the true size by a comparable factor. 



6. Angular Displacement 



In addition to changes in the observed brightness of the lensed source, the plasma lens also 
causes variations in the apparent angular position of the source. Consider an observer at a point 
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u' , who observes a ray that, as a solution to equation (|T7|), comes from the point on the lens 
plane. In the small angle approximation, the angular displacement of the source from its "true" 
(unlensed) position is {u' — Uk)6e, cf. equation (|5|). The value of {u' — Uk) depends on u and the 
parameters a and /3s- For simplicity, we will restrict our discussion to the case of a point source, 
for which /3s = 0. 

Figure ^ shows the general behavior of the angular displacement of a background source 
as viewed through a plasma lens. While this figure was computed for one particular value 
of a, the qualitative form of the angular displacement is independent of a. The figure shows 
the function {u' — Uk), for valid solutions of equation (|T7| ) as a function of u' . The lens axis 
intersects the observer's axis at u' = 0. The vertical axis, both in contour plot form and in a 
pseudo-three-dimensional form, is the relative intensity of the image coming from a particular 
value of (n' — Uk)- At values of u' far from the lens axis, e.g., u' = 5, {u' — Uk) = 0, and there is no 
angular displacement (Image A in the figure). Closer to the lens axis, the outer caustic is formed, 
e.g., near u' = 3.1. In Fig. ^, the caustic is evident by the sudden formation of a second, bright 
image at a new value of {u' — u^)- Image A also remains; that is, (u' — u^) as a function of u' 
becomes multi- valued. The new bright image is angularly displaced from the "true" position of 
the background source as evident from the non-zero value of (n' — Uk). 

Even closer towards the lens axis, the image from the caustic breaks up into two separate 
images (B and C in the figure), and their intensities diminish. Image B moves towards the original 
image (Image A) as u' decreases, i.e., as one moves closer to the lens axis. The intensity of 
Image C diminishes quickly. 

Somewhat closer to the lens axis, the inner caustic forms, e.g., u' ~ 2.6, where Images A 
and B merge into a single very bright image. Then, Images A and B disappear altogether, leaving 
only the weak Image C. This image moves towards zero angular displacement as the observer 
nears the lens axis, u' = 0. On the axis. Image C has zero angular displacement and it reaches its 
minimum flux density, 1/(1 + a) of the nominal source intensity. 

The magnitude of the angular displacement can be estimated from equation (|5|) . The angular 
displacement in radians for an image from a point on the lens plane reaching the point u' on 
the observer plane is 

56 = {u' -uk)9e, (24) 

where we have assumed that the lens axis passes through u = u' = 0. For a Gaussian lens, the ray 
path description, equation ([ITD, gives u' in terms of the solution(s) Uk, 

{u' - Uk) = Ukae-'^K (25) 

The maximum value of the difference («' — Uk) occurs for Uk = l/\/2; which is a solution to 
equation ([T7|) when the observer is located at the point u' = (ae~^/^ -|- l)/\/2. The total separation 
between the two outer caustics is twice this value, from which equation (22) is derived. This 
expression is true for all a. 
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In Fig. ^ it is evident that the maximum angular displacement occurs as the observer crosses 
the outer caustic. Images B and C are merged into one image there. Since the image is coming 
from the point Ufc = l/\/2 and the observer is at the point u' = (ae"^/^ + l)/\/2, the angular 
displacement of the merged image is 

56 = {u' -Uk){a/D), 

= -^a6ee-^/\ (26) 



7. Application to Extreme Scattering Events 

We have presented a quantitative analysis of the refractive properties of a specific optical 
system: a one-dimensional interstellar plasma lens. The specific case of a Gaussian distribution 
of electron density transverse to the line of sight has been considered. We do not believe that 
such ideal systems exist in reality, but the results of our analysis provide a starting point for a 
semi-quantitative analysis of extreme scattering events towards distant radio sources. 

Fiedler et al. (1994a) summarized the variety of light curves observed during ESEs. The 
scattering events are pronounced at lower frequencies (« 1 GHz), but usually weakly detected 
or not detected at higher frequencies (~ 10 GHz), presumably due to a combination of weaker 
ISM scattering at shorter wavelengths and intrinsic noise in the radio light curves. One exception 
to this general rule is the first detected scattering event, towards the quasar 0954+658, which 
occurred in 1981. At 8.1 GHz, the most prominent features in the light curve are four strong 
spikes, reminiscent of caustic spikes due to a point-like source being refracted by a strong lens. 

The lack of caustic spikes in most of the identified ESEs, however, implies that the typical 
plasma lenses that produce ESEs are fairly weak (small a), so that multiple imaging does not 
occur; the lenses have much smaller angular diameters than the background sources (large /3s), 
so that caustic surfaces are less pronounced; or a combination of both effects. Fey, Clegg, & 
Fiedler (1996) present multi-epoch, multi- wavelength images of sources in which an ESE has 
been identified at times when the sources were not undergoing an ESE. They find the sources 
to be compact, with typical FWHM diameters of approximately 1 mas at wavelengths between 
3.5 and 18 cm. Figure |5| indicates that, if the lack of caustic spikes is because the lenses have 
smaller angular diameters, the typical lens diameter would be 6^ ^ 0.1 mas (0.1 AU at 1 kpc). 
We therefore favor the weak-lens model as the most plausible explanation for the general lack of 
caustic spikes in ESEs. 

The best-studied ESE is that for the extragalactic source 1741—038 which occurred in 1992. 
High quality light curves were obtained at 2.25 and 8.1 GHz during the event and refraction effects 
were detected in both. Radio polarization data obtained during the event and 1.5 yr afterwards 
have been used to constrain the magnetoionic structure of the lensing medium (Clegg et al. 1996). 
VLBI observations of 1741—038, also obtained outside the ESE, have been used to determine the 
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unlensed angular structure of this source. Fey et al. (1996) used these data m conjunction with 
the Fielder et al. (1994a) model to infer that the angular diameter of the lens and the source must 
be comparable (ss 0.5 mas). 

We have chosen the ESEs towards 0954+658 and 1741—038 for detailed analysis with respect 
to our Gaussian lens model. This choice was made due to the apparent formation of caustics 
during the 0954+658 event and the quality of the data obtained during the 1741—038 event. In 
the following discussion, we will obtain semi-quantitative estimates of lens parameters and indicate 
how additional data would be useful in improving our estimates. 



7.1. 0954+658 

The 2.25 and 8.1 GHz light curves obtained for 0954+658 are shown in Fig. ^. The ESE, 
which occurred in 1981, is shown in Fig. |9| on an expanded scale. The 2.25 GHz ESE is marked 
by a flat-bottomed minimum, at which the flux density of the source is one-half its nominal 
(unlensed) value, surrounded by two peaks, which we identify below as the inner two caustics. 
The 8.1 GHz ESE is marked by a 2 Jy peak and three 1 Jy peaks bracketing an approximately 
1.5 month minimum in the source's flux density. The minimum is less pronounced at 8.1 GHz, as 
the source's flux density during the event decreases to 70% or so of its nominal value. Below, we 
identify these four, approximately 1 Jy peaks as the caustics. 

The flat-bottomed minimum in the 2.25 GHz ESE indicates that, at this frequency, the 
intrinsic angular diameter of 0954+658 is substantially smaller than the angular diameter of the 
plasma lens, cf. Fig. ^. However, the fact that the flux density during the ESE only decreases to 
one-half its nominal value can be reconciled with our lens model only if the light curve is the sum 
of a lensed and an unlensed component. If we take the amplitude of the light curve during the 
minimum to be the flux density of the unlensed component, the light curve can be represented 
by a two-component source model in which a component of flux density 0.35 Jy is lensed by the 
intervening medium and a second component of flux density 0.3 Jy is unaffected by the lens; the 
angular diameter of the lensed component also has to be sufficiently small enough to produce the 
flat-bottomed minimum. 

The unlensed component could be either a compact component of 0954+658, such as a jet 
component, whose line of sight was unaffected by the lens, or emission on larger scales that 
is completely resolved by VLBI arrays, but is unresolved to the much smaller Green Bank 
Interferometer used to obtain the light curve. Both scenarios are consistent with the observed 
structure of 0954+658. 

Gabuzda et al. (1992, 1994) conducted 5 GHz observations between 1987 and 1989, and found 
the source to consist of a compact core (< 0.5 mas) with a jet extending approximately 5 mas 
to the northwest. Within the jet they identified multiple, compact (< 1 mas) knots. With the 
possible exception of the knot KO ( |Gabuzda et al. 1994 ), however, none of the knots they identified 
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would have been present at the time of the ESE. The proper motion of the knots is approximately 
0.4 mas yr~^, large enough that most of the knots were not ejected until after the ESE (epoch 
1981.1), unless the knots have undergone substantial acceleration since their ejection. Pearson & 
Readhead (1988) included this source in a finding survey conducted in 1978. Though they did not 
image it, they found it to be compact and poorly modelled as a single gaussian component. These 
VLBI observations bracket the time of the ESE and found indications of compact, yet complex 
structure. We therefore conclude that it is plausible that, at the time of the ESE, the source 
consisted of multiple compact components, whose typical angular scales are milliarcseconds. 

On intermediate to larger angular scales, tens of milliarcseconds to arcseconds, there is 
also emission from this source. Gabuzda et al.'s (1992, 1994) VLBI observations also included 
simultaneous VLA observations. They found that approximately 15% (~ 0.15 Jy) of the flux was 
resolved out by the VLBI observations. Also, using lower resolution VLA observations Kollgaard 
et al. (1992) find a jet extending approximately 4" to the south of the VLA core; the jet's flux 
density is approximately 0.22 Jy. This extended emission would lie entirely within the synthesized 
beam of the Green Bank Interferometer, but would be unaffected by the passage of a Gaussian 
lens in front of the source. In summary, this source displays emission on a range of angular scales, 
from submilliarcsecond to arcseconds. A combination of small-scale components, whose line of 
sight the lens did not intersect, and emission on larger scales could provide the 0.35 Jy of unlensed 
flux. 

To obtain quantitative estimates of the lens' parameters, we proceed as follows. We obtain 
approximate values for a and Ps by appeal to the light curves. We refine our estimates for a and 
Ps by producing simulated light curves with slightly different values of a and /?s and comparing 
visually the simulated and actual light curves. Due to the simple model we are considering here, 
we judged a more exhaustive search of the (a, /3s) parameter space not to be worthwhile. 

Comparing the shape of the inner caustics and minimum of the 2.25 GHz ESE to Fig. ^, we 
estimate /?s to be 0.1 < /3s < 1. From the separation of the inner and outer caustics in the 8.1 GHz 
ESE and equation (^), we estimate a ~ 10 at this frequency. Since a oc A^, at 2.25 GHz we have 
a ^ 130. 

Varying the parameters and comparing simulated light curves to the actual 2.25 GHz light 
curve, we find reasonable agreement for a = 160 and /3s = 0.4. In Fig. ^ we have superposed the 
simulated light curve for this set of parameters on the 2.25 GHz light curve of 0954+658. In doing 
so, we have scaled the u' axis arbitrarily since it depends on the unknown relative transverse 
velocities of the observer and lens. The agreement with the ESE minimum is good. The simulation 
underestimates the amplitude of the inner pair of caustics. The outer caustics, for a = 160, are 
far outside the time range of the observations plotted in Fig. ^. We have compared our simulation 
to the full time range of the observations of 0954+658, Fig. |8[ We do not see any evidence for the 
existence of the outer caustics, but the predicted amplitude of the outer caustics is significantly 
smaller than the predicted amplitude for the inner caustics. The outer caustics are probably lost 
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in the noise in the hght curve. The ESE occurs during a general decrease in the source's flux 
density, probably caused by intrinsic variations. This general decrease in the source's flux density 
explains the lack of agreement between the amplitude of the actual light curve and the model light 
curve after the ESE. 

To compare our model with the 8.1 GHz light curve, we scale the light curve fit to the 
2.25 GHz ESE to 8.1 GHz. We believe such a simple scaling is justified within the context of our 
1-D simulations of a Gaussian lens passing in front of a single compact component of 0954+658. 
Below we discuss briefly the consequences of relaxing the assumptions in our simple model. 

We assume that the angular size of the lensed component of 0954+658 scales as A, as is 
generally appropriate for flat-spectrum extragalactic radio sources. Though the refractive index 
within the lens does change with wavelength, the physical size of the lens should not. Thus, the 
relative angular extent of the lensed component is (2.25/8.1) = 0.28 times smaller at 8.1 GHz as 
compared to 2.25 GHz or [3^ = 0.11 at 8.1 GHz. Our estimated value of a at 2.25 GHz, a = 160, 
corresponds to a = 12 at 8.1 GHz. The resulting model light curve is superposed on the actual 
light curve in Fig. ^. The best match for the flux densities of lensed and constant components at 
8.1 GHz is 0.15 and 0.45 Jy, respectively. 

We identify tentatively the lensed component of 0954+658 as a jet component, not the core. 
The difference in flux density at the two frequencies for the lensed component, 0.35 Jy at 2.25 GHz 
to 0.15 Jy at 8.1 GHz, indicates that the lensed component has a steep spectrum, S oc v^^-"^ . 



The emission from the source at 5 and 8 GHz is dominated by compact components (Gabuzda 



et al. 199^ , 1994; Gabuzda fc Cawthorne 1996 ; unpublished VLBI observations), generally with 



flat spectra (viz. Fig. 8). However, one of the two knots in the jet, knot K3 ( Gabuzda et al. 1994 ; 



Gabuzda fc Cawthorne 1996 ) has a similar spectral index with a 5 GHz flux density of 0.2 Jy. 
We re-iterate that, because of proper motions, none of the knots seen at 5 GHz by Gabuzda et 
al. (1994) or at 8 GHz by Gabuzda & Cawthorne (1996) would have been present at the time of 
the ESE. Nevertheless, we believe the existence of such a jet component at a later date suggests 
that such a component could have been present during the ESE. Clearly VLBI observations 
during and after an ESE would be quite useful in assessing which component(s) was lensed. 

Qualitatively, Fig. ^ shows that our simulation reproduces the four caustic spikes in the 
8.1 GHz light curve. However, the time scale for the scattering event derived from our model 
differs substantially from the observed data. The model predicts a duration of approximately 0.4 
yr, while the observed duration is about half that. Equation (22) and Fig. ^ show that the distance 



(and therefore time scale) between the inner pair of caustics is a very weak function of a in our 
model, and, therefore, a very weak function of observing wavelength. In essence, the time scale 
between the inner caustics at 2.25 and 8.1 GHz should be nearly identical. In the observational 
data, that is not the case. 



The most likely explanations for the discrepancy between our model and the data are 
due to the limitations of a one-dimensional lensing simulation. The wavelength scaling of the 
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timescale of the lensing event derived in the 1-D model does not account for the possibiUty that 
the mid-point of the lens may not cross the background source. Instead, only a chord of the 
lens, rather than the full diameter, may pass in front of the background source. Combined with 
complex wavelength-dependent structure in the background source, meaning that the predominant 
contribution to the flux density may shift in position with wavelength, the net effect will be a 
complicated wavelength scaling of the time scale for the lensing event. 

We use our estimates of a and (3s to constrain the physical properties of the lens. Since 
the better comparison is obtained at 2.25 GHz, we will use the parameter values found at that 
frequency: a = 160, /3s = 0.4. From equation ([Tsl), the maximum electron column density through 
the lens is 

A^o = 0.28pccm-3 aX-'^DO'lp-'^ . (27) 

Here is the angular size of the background source in units of milliarcseconds, and we have used 
the fact that j3s = 0s/6(^ from equation (13). 



The unknown quantities in equation ( |27D are the angular diameter of the background source, 
^s, and the distance to the lens, D. The angular diameter is ^ 1 mas ( pabuzda et al. 1992| , 1994; 



Gabuzda &: Cawthorne 1996| ; unpublished VLBI data). Toward 0954-1-658 is Galactic Loop HI, 



which Fiedler et al. (1994b) identify as responsible for the ESE. Its distance is estimated to be 
0.15 kpc. More generally, the scale height of the Galaxy's free electron layer is approximately 
1 kpc ( pTaylor fc Cordes 1995j ). At the Galactic latitude of 0954-1-658 (6 = 43? 1), the maximum 



distance to the lens would be roughly 1.5 kpc. 

With these estimates for 9^ and D, we compute a nominal column density through the lens of 
A'^o — 0.24 pc cm~^ corresponding to a distance of 0.15 kpc and an upper limit of Nq <^ 2.4 pc cm~'^ 
corresponding to a distance of 1.5 kpc. Our estimates for 9s and /3s also provide constraints on 
the size of the lens. We estimate 9e « 2.5 mas, corresponding to a 0.38 AU (3.8 AU) for a 
lens at a distance of 0.15 kpc (1.5 kpc). The free electron density within the lens is Ue ~ Nq/q 
or rie ~ 10^ cm~'^; in our model, the electron column density estimate is distance independent. 
The mass of the lens is ~ uipriea'^, where mp is the mass of the proton; for a lens at 0.15 kpc 
(1.5 kpc), Mi ~ 1.3 X 10^0 g ~ 6.5 X lO^^"^ Mq (1.3 x 10^3 g ~ 6.5 x lO^^^ Mq). 

The angular displacement of the position of the source during the ESE can be estimated from 
equation (|26|), 69 ~ 250 mas at 2.25 GHz. Phase-referenced VLBI observations can obtain absolute 
position information. Angular displacements of this magnitude should be detectable easily, even 
if the phase-referencing source is a few degrees away ( [Beasley &: Conway 1995 ). In the case of 



0954-1-658, the radio reference frame source 0951-1-693 (M81) is roughly 4° away. Of course, if the 
source is not lensed completely, relative positions between components can be determined with 
high accuracy without recourse to phase-referenced observations. We stress that the maximum 
angular displacement, and the optimal time for such observations, occurs during passage through 
the outer caustic. 



We close this section speculatively. First, we use the linear lens size and the time scale of 
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the ESE to estimate the velocity of the lens. From equation (^), v ~ 4.3a/At. The time scale 
between the inner caustics of the ESE is At ~ 0.25 yr. Thus, the velocity of the lens is 30 km s~^ 
(300 km s~^) for a lens at 0.15 kpc (1.5 kpc). Second, the lens responsible for the 0954+654 ESE 
is required to be a factor of a few larger than the source; Ps = 0.4 implies the angular diameter 
of the lens is a factor of 2.5 greater than the angular diameter of the lensed source component. 
However, the lens' diameter is also comparable with the size of the jet imaged by Gabuzda et 
al. (1992, 1994). If the lensed component was indeed a jet component, the lens would have had to 
pass in front of the jet component without affecting the fiat-spectrum core. We therefore suggest 
that the lens was anisotropic, with the major axis corresponding to the derived lens diameter. The 
minor axis would have to be smaller than the typical knot separation in the jet, probably less than 
1 mas, suggesting an axial ratio of at least 2.5:1. We simply note that anisotropic structures in 
the interstellar plasma are not unexpected (e.g., Higdon 1984, 1986). Finally, the smooth, single 
minimum in the 2.25 GHz light curve suggests that the lens encountered only a single component 
in the jet, i.e., moved northeast-southwest (or vice versa), roughly perpendicularly across the jet 
axis. 



7.2. 1741-038 

The complete light curve for 1741—038 is shown in Fig. 11^ with the ESE, which occurred in 



1992, shown in Fig. 11 on an expanded scale. At both frequencies the ESE displays a rounded 
minimum, with the 2.25 GHz ESE causing an approximately 50% decrease in the source's flux 
density and the 8.1 GHz ESE causing an approximately 30% decrease in the source's flux density. 
The lack of a flat bottom and the apparent lack of caustic spikes at either frequency suggest that 
1741—038 was only weakly lensed, i.e., no multiple images were formed and a < amin- 

Our analysis for 1741—038 proceeded in the same fashion as that for 0954+658. Our best 
model comparison to the 2.25 GHz light curve gives the following parameter values: a = 2, 
/3s = 1.0, with a single-component (lensed) source model of flux density 2 Jy. The model is 
superposed on the observational data in Fig. |l^ and replicates the minimum of the 2.25 GHz event 
very well. It overestimates the amplitude of the maxima surrounding the minimum, particularly 
the one preceding the decrease in the source's flux density, but the amplitudes of these maxima 
are not too much larger than the level of (presumably) intrinsic fluctuations in the light curve. 

Scaling the value of a to 8.1 GHz, we find a = 0.14. Scaling the source size by A gives 
/3s — 0.28, but Fey et al. (1996) show that the dominant component of 1741—038 is unresolved 
at this frequency. We therefore assume that /3s = at 8.1 GHz (in practice, there is very little 
difference between /3s = 0.28 and /3s = light curves). The single component model with a = 0.14, 



/3s = 0, and lensed source flux density of 2.55 Jy is shown in Fig. |ll|. The observed depth of 
the light curve is greater than that reproduced by the model and the observed minimum occurs 
slightly before that predicted by the model. These deficiencies may indicate that the lens did 
not pass directly over the source, crossing only at a grazing incidence, that the electron column 
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density within the lens is more comphcated that the Gaussian form we have assumed, or that the 
source contains unresolved substructure. 

Following our analysis for 0954+658, we require the source size and distance to the lens in 
order to constrain Nq, equation (p7|). Fey et al. (1996) show that the angular size of 1741—038 is 
approximately 0.5 mas at 2.32 GHz. Since /3s = 1, the lens must also be approximately 0.5 mas 
in diameter. A number of sources displaying ESEs, including 1741—038, are seen along the edges 
of radio Loop I ( [Fielder et al. 1994a ). Berkhuisjen (1973) estimates a distance of 130 it 75 pc to 



Loop I. If we assume the ESE arises from a structure associated with Loop I, the resulting peak 
electron column density through the lens is Nq ~ 10^^ pc cm^^. If we do not associate the ESE 
with Loop I, we can place a limit on its distance only by requiring that the lens be within one 
scale height of the free electron layer of the Galaxy. This distance limit is D ^ 4.5 kpc, for a peak 
column density of Nq ^ 3.6 x 10~^ pc cm~^. Although a connection between Loop I and the ESEs 
is suggestive, the line of sight to 1741—038 (and other, nearby sources which have undergone an 
ESE) passes through the inner Galaxy and some other Galactic structure may be responsible for 
the ESEs. 

The inferred size of the lens is 0.065 AU (^ 2.25 AU) if the lens is (is not) associated 
with Loop I. The density within the lens is Ue ~ 300 cm^'^, for a corresponding mass of 
1.6 X 10^^ g ~ 8 X 10"^^ M0 (6.8 x 10^^ g ~ 3.4 x 10"^^ Mq) if the lens is (is not) associated with 
Loop 1. The time scale between the inner caustics is 0.4 yr. The corresponding velocity of the lens 
is 3 km s~^ (^ llOkms"^) for a lens associated (not associated) with Loop I. 

These lens parameters are comparable to those derived for the ESE for PSR B1937+21 
(ICognard et al. 1993| ). iVo = 6.5 x 10"^ pc cm'^ and a « 0.07 AU. (They actually fit for two 



lenses, we quote the total electron column density through both lenses and the mean size of the 
two lenses.) 

As we indicated in §|6|, angular displacement occurs regardless of the value of a. Thus, 
even though caustic surfaces were not formed during lensing of 1741—038, we can still apply 
equation (pB]). We find the maximum angular displacement during the ESE to be approximately 
0.4 mas at 2.25 GHz and 0.03 mas at 8.1 GHz. It is unlikely that these angular displacements 
could be measured using phase-referenced observations. Relative position shifts of this magnitude 
between lensed and unlensed components might be detectable at 2.25 GHz, particularly if the 
source structure could be determined soon after the ESE so that proper motion of the components 
is negligible. 

Clegg et al. (1996) conducted polarization observations of 1741—038 with the VLA during 
and after the ESE in an effort to detect changes in Faraday rotation as a result of the lens passing 
across the line of sight. No rotation measure change was detected to a level of ARM ^ 10 rad m~^. 
The rotation measure change is given by 

AiiM ~ 0.81(B||)iVo, (28) 
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where is the mean component of the magnetic field within the lens parallel to the line of 
sight in microgauss. Our derived value of Nq is consistent with the lack of observed change in RM 
since the mean parallel magnetic field component would have to be (^y) ^ 120 mG (3.5 mG), if 
the lens is (is not) associated with Loop I, in order for a change in RM to have been detected. 
These magnetic field levels are 10^-10^ times larger than those typical of the interstellar medium. 
If the lens is not associated with Loop I, Clegg et al. (1988) have shown that magnetic field levels 
B ^ 2 mG can be encountered as the result of strong interstellar shocks; however, if the lens is 
associated with Loop I, the magnetic field required to produce an observable change in RM is still 
a factor of 60 larger than what can plausibly be produced. 



8. Conclusions 

We have presented the geometric optics for refraction by an interstellar plasma lens, with 
specific application to a lens with Gaussian profile of free electron column density. We have shown 
that the one-dimensional refractive properties of a lens can be characterized completely by two 
dimensionless parameters. The first parameter characterizes the refractive power of the lens and is 

where A, A'o, D, and a are the wavelength of observation, maximum free electron column density, 
distance from lens to observer, and size of the lens transverse to the line of sight, respectively. The 
second parameter characterizes the extent to which the lens' effect is diminished by the intrinsic 
size of the source and is 

Ps = 0s/Oi, (13d) 

which is the relative angular size of the background source, ^s, compared to the angular size of the 
lens as seen by the observer, 9i. 

The effect of the plasma lens is to enhance or reduce the observed brightness of the background 
source due to focusing or spreading of ray bundles, and to change the apparent position of the 
background source due to refraction. The minimum brightness of the background source occurs 
directly on the lens axis. For a point source, /3s = 0, the on-axis brightness will be 1/(1 + a) of its 
nominal (unlensed) brightness. 

On either side of the lens axis, at a distance determined by the value of a, the background 
source will have a maximum intensity due to focusing of ray bundles. If a > amin — 2.25, the 
focusing of ray bundles will become strong enough for ray paths to cross, and multiple imaging 
of the background source will occur. For the case of a Gaussian lens profile, an observer located 
within a region of multiple imaging can see as many as three images of the background source, 
each with a different brightness and appearing to come from slightly different directions. The total 
brightness of the source, obtained by adding the brightness of the multiple images, will be greater 
than the nominal brightness of the source. 
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Separating the multiple imaging regions from the single image regions are caustic surfaces, on 
which two of the three multiple images blend together and, in the limit of geometrical optics, the 
observed brightness of the background source grows without bound. An observer passing from a 
single-image region through the multiple imaging regions and back to a single-image region can 
use the observed position of the caustics to determine the magnitude of the parameter a. The 
separation between the inner caustics can also be used to constrain the diameter of the lens. 

The primary effect of increasing the angular size of the background source in relation to the 
angular size of the lens, i.e., increasing (3s, is to smooth out the light curve and reduce the observed 
amplitudes of maximum and minimum intensity. 

We have generated sample light curves from our model and compared them to extreme 



scattering events observed towards the extragalactic sources 0954-1-658, §7.1, and 1741—038 



§7.2. These two sources were chosen because their ESEs have the highest quality data available 
and, in the case of 1741—038, polarization observations also exist. In general we find reasonable 
agreement between the observed and modelled light curves at 2.25 GHz. We have far more 
difficulty recreating the 8.1 GHz light curves by utilizing only the appropriate wavelength scalings 
for a and /5s. The discrepancies between the modelled and observed light curves may result from 
some combination of substructure within the lens, an anisotropic lens shape, a lens which only 
grazes the source rather than passing completely over it, or unresolved substructure within the 
extragalactic sources. 

For 0954-1-658 we find it to be strongly lensed at 2.25 GHz. Interpreting the light curve for 
this source is difficult because of uncertainties about the source's structure at the time of the ESE. 
We favor an interpretation in which the lens passed over a component in the source's jet. The 
lens toward this source is probably associated with radio Loop HI so that the lens was 0.15 kpc 
distant. The peak column density in the lens was 0.24 pc cm~'^, its size was 0.38 AU, the electron 
density within the lens is 10^ cm~^, and its mass was 6.5 x 10~^^ ^0- The lens also caused the 
source's position to wander by as much as 250 mas at 2.25 GHz, an easily detectable amount 
using phase-referenced VLBI observations or by measuring the relative positions between source 
components. 

For 1741—038, the ESE was caused by a weak lens. The lens toward this source is likely to 
be associated with radio Loop I and at a distance of 0.13 kpc. The lens' peak column density 
is 10~^ pc cm~^, size is 0.065 AU, density is 300 cm~'^, and mass is 10"^^ Mq. From Clegg et 
al.'s (1996) upper limit on the change in Faraday rotation towards 1741—038 during the ESE, we 
place an upper limit of approximately 100 mG on the magnetic field within the lens. The angular 
position wander caused by this lens was only 0.4 mas at 2.25 GHz. An angular displacement 
of this magnitude would not have been detectable for 1741—038; many other sources that have 
undergone ESEs have more complicated structures than 1741—038, however, and relative position 
shifts of this magnitude between lensed and unlensed components might be detectable in those 
sources. These lens parameters are comparable, within 50%, to those derived for the ESE for 
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PSR B1937+21 ( ICognard et al. 1993| ). Nq = 6.5 x 10~^ pc cm-^ and a ^ 0.07 AU. (They actually 



fit for two lenses, we quote the total electron column density through both lenses and the mean 
size of the two lenses.) 

A key assumption of our plasma lens explanation for ESEs is that the lenses are discrete 
objects. The lens properties we have derived are consistent with this assumption. If we 
assume the plasma temperature within the lenses to be 10^ K, the pressure within the lenses is 
nT ~ 3 X 10^-10^ K cm~^. These pressures are well in excess of the average ISM pressure of 



roughly 4000 K cm ( Kulkarni fc Heiles 1988 ). Such lenses would be either highly transitory 



features or embedded in high-pressure environments. Romani et al. (1987) suggest ionization 
fronts, cooling instabilities, or both associated with old supernova remnants as sites of lenses. 

We emphasize that the difficulties in interpreting the light curve for future ESEs could be 
ameliorated considerably by VLBI imaging of a source as soon as possible after, or indeed during, 
an ESE. Considerable angular displacements might also be detected with such observations. 
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Fig. 1. — Example of ray path geometry through a plasma lens, (a) The ray from an infinitely 
distant point source is incident normal to the lens plane, and strikes the lens plane at the coordinate 
X. Because of the non-uniform free electron column density Nc{x), the lens refracts the ray through 
the angle 9^- The ray travels the distance D between the lens plane and the observer's plane, and 
strikes the observer's axis at the point x' . (b) The same situation occurs, except that the ray is 
incident with an oblique angle 6i upon the lens plane, i.e., the infinitely distant background point 
source is off- axis. 



Fig. 2. — Schematic diagram of refraction by a Gaussian plasma lens. See §4.1 for a complete 
description of this figure. 

Fig. 3. — Light curve of an infinitely distant point source refracted by a lens with a = 36. 

Fig. 4. — Light curves produced by refraction through lenses with a range of values of a, as defined 



by equation (15). 



Fig. 5. — Light curves produced by an a = 25 lens, with various relative angular sizes of background 
source and lens. The parameter /3s is the angular size of the (unlensed) background source in units 
of the angular size of the lens, as seen by the observer. 

Fig. 6. — Separation between the outer caustics, Au^; between the inner caustics, Au[; and between 
the pairs of caustics on either side of lens axis, (Au^ — Au[)/2, as functions of the parameter a. 

Fig. 7. — Brightness and angular displacement of images formed by the Gaussian plasma lens as 
functions of observer position u'. See §|6| for a complete description of this figure. 

Fig. 8. — The complete light curve of the quasar 0954+658 during the period of observation from 
1979-1994. The scattering event occurs at approximately 1981.1. 

Fig. 9. — The light curve of the quasar 0954+658 during an extreme scattering event. Panel (a) 
shows the data obtained at 2.25 GHz, compared to our lens model (dashed line) with a = 160 and 
Ps = 0.4. A two-component model for 0954+658 was used: a lensed component of fiux density 
0.35 Jy, and an unlensed component of fiux density 0.3 Jy. Panel (6) is the 8.1 GHz data compared 
to a two component model with a = 12, /3s = 0.11, lensed flux density 0.15 Jy, and unlensed flux 
density 0.45 Jy. 

Fig. 10. — The complete light curve of the quasar 1741—038 during the period of observation from 
1983-1994. The scattering event occurs at approximately 1992.4. 
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Fig. 11. — The light curves of 1741—038 during an extreme scattering event, (a) The dashed hue 
is our model comparison with a = 2, f3s = 1, and a flux density for 1741—038 of 2 Jy at 2.25 GHz. 
(6) The model uses a = 0.14, = 0, and assumes an 8.1 GHz flux density of 2.55 Jy. 
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